Einstein-Podolsky-Rosen entanglement in two-well BEC systems: ground states 
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We consider how to generate and detect Einstein-Podolsky-Rosen (EPR) entanglement and the 
EPR-steering paradox between groups of atoms in two separated potential wells in a Bose-Einstein 
condensate (BEC). We present experimental criteria for this form of entanglement, and propose 
experimental strategies using adiabatic cooling to the ground state. These approaches use either 
two or four spatial and/or internal modes. We also present higher order criteria that act as signatures 
to detect the multiparticle entanglement present in this system. We point out the difference between 
spatial entanglement using separated detectors, and other types of entanglement that do not require 
spatial separation. The four-mode approach with two spatial and two internal modes results in an 
entanglement signature with spatially separated detectors, conceptually similar to the original EPR 
paradox. 



I. INTRODUCTION 

The Einstein-Podolsky-Rosen (EPR) paradox [l[ es- 
tablished a link between entanglement and nonlocality 
0| in quantum mechanics. The extent to which entan- 
glement can exist in spatially separated macroscopic and 
massive systems is still essentially unknown. Entangle- 
ment in optics however has been extensively studied and 
numerous experiments have shown evidence for it 043 ■ 
An important distinction is that optical entanglement in- 
volves (nearly) massless particles, and hence is a much 
less rigorous test of any gravitational effects present. 

Generation of EPR entanglement between two mas- 
sive systems therefore represents an important challenge. 
Such entanglement is a step in the direction of funda- 
mental tests of quantum mechanics, and is relevant to 
the long term quest for understanding the relationship 
between quantum theory and gravity. Ultimately, one 
would like to demonstrate spatially entangled mass distri- 
butions, and this appears much more promising for ultra- 
cold atoms than for room-temperature atoms. For this 
reason, we focus on ultra-cold BEC environments here. 
This is also relevant if BEC interferometry is to be use- 
ful to those areas of quantum information and metro logy 
where entanglement is known to give an advantage [8l-ll3| . 
In this paper, we study strategies for generation of EPR 
entanglement between Bose-Einstein condensates (BEC) 
confined to two spatially separated potential wells. 

Quantum correlations and EPR tests for Bose-Einstein 
condensates have been suggested previously, with strate- 
gies involving molecular down-conversion [I4| | and four 
wave mixing interactions [T§ - fl7l |. among others. Early 
experiments measuring free-space correlations demon- 
strated promising signatures of increased fluctuations as- 
sociated with entanglement , but were unable to 
conclusively demonstrate entanglement or squeezing via 
reduced fluctuations, largely due to measurement ineffi- 



* pdrummondOswin. edu.au 
t mdreid@swin.edu.au 



ciencies. This has improved with recent multi-channel 
plate detection methods, but detection efficiency still re- 
mains an issue [20]. Entanglement has also been mea- 
sured, very recently, for distinct but nearly spatially su- 
perimposed modes |2lM23| in an optical lattice. 

Here, we are motivated to study the two well case, 
in view of experiments that have used this or similar 
systems to confirm both sub-shot noise quantum cor- 
relations [24|], and multiparticle entanglement among a 
small group of atoms (25l . [26| . For much larger num- 
bers of atoms (~ 40, 000), nearly quantum limited inter- 
ferometry has been recently verified J27|, showing that 
trapped atom interferometry has the potential to reach 
mesoscopic sizes. There have also been a number of pre- 
vious theoretical studies [2^, [29| that outline different 
proposals and entanglement signatures. 

The goal of this paper is to clarify what it means to 
have an EPR entanglement between groups of atoms in a 
BEC, and outline some different strategies for achieving 
this goal. We define EPR entanglement as being that 
entanglement existing, at least in principle, between two 
spatially separated systems, so that it could potentially 
realize an EPR paradox. For EPR entan glem ent to be 
claimed, three properties must be evident |30(: 

1. Two systems must be shown entangled through lo- 
cal measurements at spatially distinct locations. 

2. The nature of the entanglement criterion should be 
to confirm an EPR paradox. This requires mea- 
surement of sufficiently strong correlation between 
the two systems, for two non-commuting"EPR" ob- 
servables like position and momentum [7|. A more 
sophisticated approach would allow other entan- 
glement measures, such as those for "EPR steer- 
ing" (3ll - [36| which reveal an inconsistency between 
EPR's local realism and the completeness of quan- 
tum mechanics for other types of measurement 
strategies. 

3. To justify EPR's no "spooky action-at-a-distance" 
assumption the measurement events should be 
causally separated. 
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For large groups of atoms, the task of detecting EPR 
entanglement is much more feasible when the emphasis 
is on the EPR paradox itself, rather than on the failure 
of Bell's local hidden variable model 0- For spatially 
separated systems, the detection of sufficient correlation 
of locally defined EPR observables so that entanglement 
is confirmed [37l - [39l | would represent an achievable first 
benchmark. This by itself is not direct evidence for the 
EPR paradox, although it is a necessary condition. The 
second step of confirming the paradox has been carried 
out for photons [30j . and also appears achievable for 
atoms. The last step is probably the most difficult for 
atoms. It would require either very fast measurements 
in one vacuum chamber, or hybrid techniques involving 
two separated BECs with coupling via atom-photon in- 
terfaces |4(| , in order to achieve causally separated mea- 
surement. 

There are many possible strategies for generation of 
spatial EPR entanglement. Early experiments employed 
two photon cascades and, later, optical parametric down 
conversion, to generate entangled photon pairs 043- 
Continuous variable EPR entanglement between two 
fields, in a so-called "two-mode squeezed state" (4l| . 
was also generated using parametric down conversion 
d, 0, E3 ■ Such entanglement gave evidence for an EPR 
paradox [30j , although true causal separation of measure- 
ment events was not demonstrated in these experiments. 

The paper is arranged as follows. In Section [TT] we 
give a general introduction to the different possible en- 
tanglement strategies. Section IIIII focuses on signatures 
for demonstrating entanglement, as well as the hierarchy 
of nonlocality measures including EPR-steering and Bell 
indicators. Section IIVI considers entanglement prepara- 
tion in a two-well system, modeled as two modes with 
boson operators a and b. In this case, the S-wave 
scattering intra-well interactions, given by Hamiltonians 
H = ga) 2 a 2 and H = gb' 2 b 2 , provide the local non- 
linearity at each well, while the coupling or tunneling 
inter-well term, modeled as H = n(a^b + atf) generates 
inter-well entanglement. Here the intra- and inter-well 
interactions act simultaneously, in the adiabatic forma- 
tion of the ground state. Section [V] treats a four-mode 
generalization of this. This has the advantage that EPR- 
entanglement can be measured using atom counting at 
each site, without the use of a local oscillator. Our con- 
clusions are summarized in Section I VII with technical 
details given in the Appendices. This paper is based 
on the preliminary ideas presented in a Letter [29j . A 
second class of entanglement strategies using dynamical 
techniques will be analyzed in a subsequentpaper. 



II. ENTANGLEMENT STRATEGIES 

A. Prototype states for two-mode entanglement 

Suppose two spatially separated systems are describ- 
able as distinct modes, represented by boson operators 



a and b. There are two prototype states that one can 
consider, that can give multiparticle EPR entanglement. 
The first, which we call particle-pair generation, is cur- 
rently the most widely known and used @- We consider 
an entangled state with number correlations: 

oo 

k/>)i/ = X! c »«l n ><" ( IL1 ) 

This type of two-mode squeezed state gives two-particle 
correlations arising from a pair production process H = 
Ka^b^ + K*ab where (ab^) = but (ab) ^ 0, and the num- 
ber difference is always squeezed (43|, [44|. These EPR 
states are formed in optics with parametric down con- 
version I30B , and similarly in nondegenerate four wave 
mixing [45j. Since they are not number-conserving, they 
are not typical of states formed in coupled two-well ex- 
periments, although they have been generated in recent 
BEC experiments using spin or mode-changing collisions 
[2LM23l | . We note that a full analysis of such a system 
requires consideration of the dynamics of the coupling 
to the pump that supplies the pairs of atoms. Gener- 
ally, there is a back-reaction on the pump mode, causing 
excess phase-noise that limits amount of squeezing or en- 
tanglement that can be generated (4(| |47| . 

In this paper, we will focus on a second form of EPR 
entanglement, which we call number conserving. This 
occurs, for example, when fixed number states are input 
into a beam splitter: H = k(<2^6 + atf), so that (ab) = 
but (atf) ^ 0. We consider an entangled number- 
conserving state of form (48l - [53j : 

N 

\Tp)i = J2c n \n)a\N-n) b . (11.2) 
n=0 

This is similar to the states generated in optical fiber 
pulsed squeezing (54|, and is the closest to the state pre- 
pared in some recent two-well BEC experiments, where 
the total number is conserved 0, HH . We will examine 
situations with number anti-correlated states prepared 
through evaporative cooling and adiabatic manipulations 
near the ground state of a coupled-well system, as shown 
in Fig. □ 

B. Experimental strategies 

Before examining detailed solutions for an interacting 
BEC, it is useful to summarize how two-mode number- 
conserving entanglement can be generated, in schematic 
form. We consider how to generate entanglement be- 
tween two groups of atoms in separated potential wells in 
a BEC. What is useful is a combination of nonlinear lo- 
cal interactions to generate a nonclassical squeezed state 
in each well - together with a nonlocal linear interac- 
tion to produce the entanglement between two spatially 
distinct locations. In the case of the BEC, the S-wave 
scattering can provide a nonlinear local interaction, and 




Figure 1. (Color online) A double- well, one spin orientation 
BEC. a,b are operators for two modes at A and B. The a and 
b are prepared with a two mode number difference squeez- 
ing and an entanglement, by adiabatic cooling to the ground 
state. In Section V, we develop signatures to detect the inter- 
well entanglement, using inter-well spin operators. 
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Figure 2. (Color online) A double-well, two spin orientation 
BEC. We suppose the modes and hi are spatially separated. 
Modes oi, a 2 could be different spatial modes, or different 
spin components of the same well. Pairs a\, b\ (and ai, 62) 
can become entangled, due to the inter-well couplings. We 
allow for the asymmetric case where pair 0,2 and 62 have much 
greater numbers than a\ and 61 (N2 3> Ni) and also consider 
a case where modes 02 and &2 need not be entangled (K2 — 0). 



quantum diffusion across a potential barrier acts like a 
beam-splitter to provide the final nonlocal linear interac- 
tion. Both effects occur at the same time in the schemes 
treated here, in Sections IIVI and IVl 

We will show in Section ITVl that the entanglement gen- 
erated for the two-well ground state with a fixed num- 
ber of atoms can translate to an EPR steering type of 
entanglement [HI, [3(| (Fig. Q]). For an actual demon- 
stration of this sort of EPR entanglement, however, one 
must use signatures that involve local measurements, for 
two spatially separated observers (often called Alice and 
Bob), at sites A and B. One can use local oscillator (LO) 
measurements at each site, that provide phase shifts or 
their equivalent between the measured and LO modes 
[l6| . Here, we propose an alternative though similar four- 
mode strategy, as shown in Fig. [2] We consider two types 
of gedanken-experiment as follows: 

• Two-mode entanglement preparation then 
analysis: suppose the initial state is entangled 
through state preparation in a double-well poten- 
tial. Here the nonlinearity occurs during state 
preparation, and the other steps are to enable the 
entanglement to be transformed into an easily mea- 
surable form. Experimentally, this appears rela- 
tively simple, although obtaining measurable sig- 
natures can be a challenge if strictly local measure- 
ments are used, as in an EPR paradox demonstra- 
tion. This strategy is discussed in Section IIVI 

• Four-mode entanglement preparation then 
analysis: we consider four-mode states created 
through evaporative cooling and adiabatic process- 
ing in a double-well potential with two spin states 
in each well. Experimentally, this is more complex, 
but obtaining measurable signatures only requires 
local measurements in each well. This strategy is 
discussed in Section IVl 



In both two and four mode cases, the basic idea is: 

1. Correlated ground state preparation, through evap- 
orative cooling in a potential well with linear cou- 
pling between wells. 

2. Local Rabi rotation to a superposition of internal 
spins, thus choosing an EPR measurement angle. 

3. Measurement, usually from absorption imaging, 
giving occupation numbers. 

III. MEASUREMENT STRATEGIES 

A. Particle-pair EPR entanglement 

In the EPR proposal PJ, the paradox arose from cor- 
relations between the positions and momenta of two 
particles emitted from the same source. With optical 
or atomic Bose fields, one can define the quadrature 
phase amplitudes of the modes, as Xa = a' + a, and 
Ya = {a' — o,)/i, and similarly for mode b. These have 
the same commutators as position and momentum in a 
particle system. These quadratures, which are measured 
using local oscillators, can be used to detect the entan- 
glement of the pair-correlated states discussed above, as 
performed recently for atoms by Gross et al [22| . 

1. Entanglement 

The original EPR paradox focused on states that 
showed correlation and anti-correlation for position and 
momentum respectively. Duan et al and Simon [13, HH 
showed that entanglement between modes a and b is con- 
firmed if 

A 2 (Xa-X b ) + A 2 (P a + Pb) <4. (III.1) 
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In the case pil.l|) . the 4 arises from the commutation 
relation [a, a'] = 1, and reflects the quantum noise asso- 
ciated with the four observables. 

A full analysis of a local oscillator measurement shows 
that it is actually equivalent to a spin measurement, and 
closely related entanglement measures have been devel- 
oped for use with spin measurements. In particular, one 
can show entanglement using the spin version of (|III.1I) 

M, 



B. Hillery and Zubairy type entanglement 

The signatures given above are not the only possible 
ones, and in fact there are many signatures both for en- 
tanglement and the EPR paradox. A particularly useful 
alternative form was proposed by Hillery and Zubairy 
(56| , and Cavalcanti et al [36l [57J , as described and gen- 
eralized below. 



tJ%) + &VI±Jb)< 1(^)1 + 1(^)1, (ni.2) 

and also the (more sensitive) Heisenberg-product entan- 
glement criterion [13 



A 2 (J «+-/2 



e+Tr/2 
B 



)< 



(Jl) 



(III.3) 



2. EPR paradox entanglement 



EPR's argument assumes that one observer (Alice) can 
make precise predictions for the outcome of position or 
momentum measurements made by a second, distant ob- 
server (Bob). To quantitatively demonstrate the EPR 
paradox, the level of correlation in Alice's predictions is 
compared with the quantum limit for a local state that 
might predetermine Bob's statistics. Thus, for the EPR 
paradox, the relevant quantum noise level is that of one 
observer, B, alone, which implies a stricter condition for- 
mulated in terms of conditional variances Q. EPR para- 
dox entanglement is observed when 



A{X B \X A )A{P B \P A ) < 1 



(III.4) 



Since the precise choice of measurement at mode A is not 
important, only the inference, this criterion is sometimes 
more generally written as Ai n fX B Ai n fP B < 1. In the 
case of pil.l[) . the stricter condition 



A 2 {X A - X B ) + A 2 (P A + P B )<2 



(III.5) 



if satisfied would imply an EPR paradox entanglement 

MM- 

While the signatures (|III.HIIIL5j) have been used ex- 
tensively in relation to EPR- type entanglement, we do 
not use them in this paper. Evaluating, we find in fact 

D = A 2 (X A - X B ) + A 2 (P A + P B ) 

= 4(1 + (at, a) + {b\b)-{a, b) - (a\tf)) 

(where here we use the conventional notation that 
(x,y) = (xy) — (x)(y)). Since detection of EPR entan- 
glement using this criterion requires at least D < 4, we 
see that if (a) = (b) = 0, then (ab) ^ is a minimal 
requirement to use this signature. Where the number N 
of atoms is fixed, the ground state solution is of the form 
(III.2[) . for which (ab) — 0. For this reason, we focus on a 
different set of EPR entanglement criteria. 



1. HZ entanglement criterion 

Modes a and b are entangled if [HJ 

\(a m (6t)")| 2 > ({o)) m a m (tf) n b n ). 



(111.6) 



This criterion is useful for states with a large phase mo- 
ment (atf), as will be studied in Section V. In Ref |29| 
it was suggested how to rewrite this criterion for m = n. 
For any nonhermitian operator Z, we consider the gen- 
eralized variance, which must be nonnegative: 

A 2 Z = ((tf - (Z 1 ")) (Z-(Z))) = (Z*Z)-(Z*){Z) > 0. 

(III.7) 

Defining Z = a m (&*)" , we find it is always true (for any 
state) that 

|(a m 6 tm )| 2 - (a) m a m tf m b m ) < (a tm a m ([6 m 6 tm , 6 tm 6 m ])>. 

(III.8) 



Thus, the HZ criterion (|III.6P confirms entanglement if: 
(a) m a m tf m b m ) - |(a m 6 tm )| 2 



< Eft* = 1 



(al m a m (b m tf m - tf m b m )) 



< 1. 



(III.9) 



It is also possible to derive a criterion using the commu- 
tators for mode a. Hence the HZ entanglement criterion 
pil.9[) is best written with the optimal choice of denomi- 
nator, corresponding to the minimum of (a^ m a m (b m b' m — 
tf m b m )} or (tf m b m (a m at m - a) m a m )). 

The first order (m = n = 1) HZ criterion for entangle- 
ment becomes 



< E\ 



,,,, = ( „t<y -iff)i' 



The second order WL entanglement criterion is ob- 
tained by using the power m = 2 with the identity 
[b 2 ^ 2 , W 2 b 2 ] — 46 t 6 + 2. Entanglement is then observed 
if 



° : Ehz 1 + (at2 a 2(46t6 + 2)) < l - 



(III.ll) 

We show in Appendix B that the higher order criteria 
with m > 1 are useful for detecting a multi-particle en- 
tanglement. 
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2. HZ EPR-steering 



Criteria similar to IIII.6I that reveal not just entangle- 
ment, but the EPR-steering paradox P, 0, HU, [35| or vi- 
olation of local hidden variable theories (Bell's theorem) 
0|, have been derived in Refs [HJ and (57). In this paper, 
we consider two sites only, and focus on the entanglement 
and EPR-steering cases, since the Bell inequality derived 
in Ref [13] requires at least three sites. 

The terminology "EPR-steering" is used to describe a 
more general EPR paradox situation that can involve 
any type of measurement, not just conjugate spins or 
position/ momentum [32l l35l . The EPR paradox was 
discussed by Schrodinger [31], who pointed to the phe- 
nomenon of "steering". Criteria for "steering" can be de- 
veloped using the asymmetric local hidden state separa- 
ble model developed in Ref. J32|. Analysis tells us that 
violation of this model reveals inconsistency of local re- 
alism with the completeness of quantum mechanics, and 
is thus also a criterion for a generalized EPR paradox 

us, mum. 

An EPR-steering nonlocality is detected if 



|(o m 6 tn )| 2 > (a tm a m ( 



h \n h n + b ntfn 

2 



)). (111.12) 



The proof follows from straightforward application of 
methods is given in [36j which derived the EPR steering 
result for m = n = 1. Similar higher moment extensions 
of the CFRD Bell nonlocality [57J have been considered 
in [58]. This criterion can also be rewritten in terms of 
the HZ entanglement parameter ([III. 10|) . so that EPR- 
steering entanglement is confirmed if: 



< E 



(m) 
HZ 



= 1 



(a) m a m tf m b m ) - \(a m tf m )\ 2 1 
(a,1 m a m (b m tf m - tf m b m )) < 2' 



(111.13) 



IV. GENERATION OF NUMBER-CONSERVING 
TWO-MODE ENTANGLEMENT 

We next turn to physical means to generate and mea- 
sure entanglement and EPR-steering in two-mode phys- 
ical systems. We focus here on the gedanken-expeviment 
of Fig[U with explicit spatial separation of the two modes. 
Of course, one can use these criteria to discuss entangle- 
ment of two spin degrees of freedom as well, but unless 
these are localized with a spatial separation, there is no 
causal separation possible: and hence no EPR interpre- 
tation. 



Hillery and Zubairy [56| have written the first order 
criterion (jlll.lOj) in terms of the variances of inter-well 
Schwinger spins, defined as: 

J'Ib = (a% + atf) /2 
- (ah-atf) /(2i) 



J 



¥ 

AB 

J AB 
J 2 

J AB 

N AB 



(a) a 



rfb) /2 
iV(iV + 2)/4 
a)a + tfb. 



(IV.l) 



Where the outcomes for Nab are fixed, the spin is fixed 
as J = N/2. The resulting ground-state HZ entangle- 
ment criterion is given by Eq. (|III.10[) which can also be 
written as: 



0<E HZ 



(AJf s ) 



(AJJ B ) 



{N AB ) /2 



< 1. (IV.2) 



It should be noted here that this type of spin-operator 
variance has been measured experimentally [2~H by ob- 
serving the interference between the two modes, on ex- 
panding the atomic clouds after turning the traps off. 
However, as we discuss later, this strategy cannot be 
readily interpreted in the EPR sense, due to the lack 
of separation during measurement. 

The best entanglement as measured by pV.2[) is given 
when the sum of the two variances of Jab an< ^ ^ab ^ s 
minimized. This sum can never be zero, meaning that 
ideal entanglement of Ehz = cannot be reached, be- 
cause the spins Jab anc ^ ^ab ^° n °t commute. However, 
the sum becomes asymptotically small for large N, in 
which case large noise appears in the third spin Jab- 
The lower bound for the sum of the two variances has 
been obtained by [59]: 



(AJ 



ab) 



J 



>Cj/J 



(IV.3) 



where the coefficients Cj are given in that reference. The 
reduction of the sum (AJ^ B ) +(A Jab) below the stan- 
dard quantum limit (given by J = (Nab)/2) is referred 
to as "planar squeezing", and represents the onset of HZ 
entanglement. 

The inequality pV.3[) leads to a useful result that will 
be used later in Section HVl Since a large spin J can only 
be obtained where the number of atoms N is large, very 
small squeezing necessarily implies an entangled state 
with a large mean (N). A full analysis is given in the Ap- 
pendix A. This type of approach has been used by Soren- 
son and Molmer [60j, where measurements are made of 
the variance of Jz, and has proved useful in deducing 
multiparticle entanglement (25l . |26| . 



A. Spin-operator HZ criterion 



B. Beam splitter with fixed number input states 



It is convenient to quantify entanglement in the 
number-conserving case using spin-operator methods. 



Possibly the simplest number-conserving entangled 
state is obtained with a number-squeezed input, together 



(i 



with a beam splitter interaction 

H = K.a f b + K *ab\ (IV.4) 

which models the exchange of atoms that can take place 
between wells. 




Figure 3. (Color online) A Fock number state \N) incident 
on a beam splitter produces an entangled state ()II.2[) . 

On defining output (a, 6), input (aj„) and vacuum (a v ) 
input modes (Fig. [3J [3]) , one can write the beam splitter 
transformation as 

a = (a m + a v )/V2 (IV.5) 
b = (a in - a v )/V2. 

1. Single number-state input 

We first consider the simplest case of N atoms input 
to one port of the beam splitter (Fig.J3J ■ This is equiva- 
lent to the linear interferometer case [25| in which a fixed 
number of atoms are initially in one BEC well. These are 
then redistributed between wells via a number conserv- 
ing mechanism. Using IIV.51 the final state is number- 
conserving pi.2[) : 

N 

\out) = J2 c n\n)a\N-n) b , (IV.6) 

n=0 

where c„ = yJW\/ v / 2 N n\(N - njl. This state pV^j) 
is entangled for all N. The entanglement can be de- 
tected using the Hillery and Zubairy entanglement mea- 
sure (|III.9|) . Higher order (up to iV-th) entanglement is 
also evident using the E^\ criteria (Fig. [4}. 

This method generates a relatively small degree of en- 
tanglement, however, (Fig. [4]), and will later be com- 
pared with the much more significant entanglement ob- 
tainable using nonlinear BEC interactions. 



2, Double number state input 

We next consider a double Fock number state |AT)|iV) 
incident on a beam splitter (Fig. 0, as a model for 





\ 

. \ 

\ 

V 

* s 

- ( • S. 
- \ ^ ^ ■ 


Entanglement 


b HZ 

L HZ 

r.<3) 
£ HZ 

;^hz 


EPR steering 



10 20 30 40 
N 



Figure 4. (Color online) A single Fock number state with 
beamsplitter is entangled {Ehz < 1) by the HZ entangle- 
ment criteria, Eq. (|III.9[) . Higher order entanglement is indi- 
cated by the dashed lines. The correlation does not confirm 
EPR-steering entanglement from Eq. 13f) . which requires 
Ehz < 0.5. 




|2n^ |2(N-n)> b 

Figure 5. (Color online) A double Fock number state inci- 
dent on a beam splitter also produces a number-conserving 
entangled state. 

the case where there is initially a fixed, equal number 
of atoms in each well. 

The output state after an exchange between the wells 

is 

N 

\out) = J2 c n \2n) a \2(N - n)) b (IV.7) 

n=0 

where c„ = {-l) N - n ^f(2n)\^ (2{N - n))\/ [2 N n\(N - n)\] ■ 
In this case, entanglement is again present for all N, 
but cannot be detected via the first order entanglement 
criterion pil.iop . 

Entanglement can however be detected via the second 
order HZ entanglement criterion Eq. (jlll.llj) , which indi- 
cates an entanglement involving superposition are num- 
ber states different by two particles. The fourth-order 
entanglement E^ is also evident, indicating superposi- 
tions involving states separated by four particles. The en- 
tanglement measure E^ is sufficiently strong that EPR 
steering can also be confirmed via Eq. (|III.13j> with 
m = n = 2, as shown in Fig. [5J though this is diminished 
at higher N. 
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1.5r 



Figure 6. (Color online) HZ entanglement criterion using 
a double Fock number state and beam splitter. The graph 
shows the criterion (|III.9|) for m = n = 2 (solid blue line), 
and m — n = 4 (red dashed line). EPR-steering is observable 
with m — n = 2 and N < 5. 



C. Nonlinear case: BEC ground state 

We now examine how to enhance the entanglement 
over the linear case above, by using a local number- 
conserving nonlinearity. 



N= 100 



1.0 



N 
X 



a, b 



LU 



0.5 




Entanglement 



-2 2 4 6 8 10 12 14 16 

N g/K 

Figure 7. (Color online) Entanglement in the ground state 
of the BEC Hamiltonian Eq. (|IV.8|) . using the HZ criterion 
Eq. (|III,10[) . plotted against the coupling constant for both 
positive and negative couplings, for TV = 100 atoms. Plots 
show HZ entanglement as a function of Ng/n, with k > 
held fixed and g varied. The mean spin is in the direction de- 
fined by Jab- Ehz < 1 indicates a two-mode entanglement; 
Ehz < 0.5 indicates EPR steering. The dashed red line gives 
the HZ criterion E' HZ for the rotated modes a' and b' . The 
predictions for the respective second order entanglement cri- 
terion E%\ (Eq. are given by the dotted and starred 



1. Two-mode BEC Hamiltonian 

We solve for the ground state of a two-component BEC 
(Fig. [T|), as modeled by the following two- mode Hamil- 
tonian: 

H = K(a t 6 + a6 t ) + |[aVaa + tftfbb]. (IV.8) 

Here k denotes the conversion rate between the two com- 
ponents, denoted by the mode operators a and b. and 
g cx a^o is the nonlinear self interaction coefficient [4j|, 
proportional to the three-dimensional S-wave scattering 
length, a,3D- The first term proportional to n describes 
an exchange of particles between the two wells (modes) in 
which total number is conserved. This term is the linear 
term equivalent to that for a beam splitter. The second 
nonlinear term can be thought of as creating squeezing. 
It models nonlinearity in many systems, including disper- 
sive optical bistability, from which squeezing is known 
to result. The Hamiltonian is a general one, and this 
model also applies to other systems such as optical cavity 
modes or superconducting wave-guides with a nonlinear 
medium. 

The ground state obtained using standard matrix tech- 
niques, and only depends on the dimensionless ratio g/n. 
We consider a total of N atoms: the number in well a 
is N a — a^a and similarly, Nb = b^b. Entanglement be- 
tween the modes a and b, and hence between the two 
wells, can be detected via the first order HZ entangle- 
ment criterion Eq. 10|1 . for both attractive (g < 0) 
and repulsive {g > 0) regimes, as shown in Figs. [7] and [5] 
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Figure 8. (Color online) Same as Fig. [7] but for much lower 
particle number, with N = 6. First order entanglement in 
(a, b) is shown by the solid blue line, with second order en- 
tanglement shown by the purple line with dots. First order 
entanglement in (a', b') is shown by the dashed red line, with 
second order entanglement shown by the black line with stars. 
The second order entanglement criterion becomes more sen- 
sitive where the nonlinearity is higher in both cases. 



2. Attractive interactions 

The best entanglement (i.e. the smallest possible value 
for Ehz) is given when the sum of the two variances of 
J^b an d Jab i s minimized. As explained above in Sec- 
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tion HVl this sum can never be zero, because the spins 
do not commute, meaning that ideal entanglement of 
Ehz — cannot be reached. The lower bound for the 
sum of the two variances has been obtained by (59j and 
is given by equation pV.3p . The reduction of the sum 

(A J%b) + (^^ab) below the standard quantum limit 
(given by J = (Nab) /2) is referred to as "planar squeez- 
ing", and represents the onset of HZ entanglement. 

The absolute lower bound for Ehz is predicted for 
the BEC ground state of (|IV.8j) in the attractive regime 
where g is negative, as shown for N = 100 in Fig. 
and for N = 6 in Fig. [p This critical case has been 
studied and explained in |59l) and We note however 
that the minimum Ehz becomes asymptotically small for 
large N, in which case large noise appears in the third 
spin Jab- The degree of entanglement increases with N 
the number of atoms, according to the relation for Cj 
obtained in (59| . 

For a given number of atoms N, the best possible HZ 
entanglement is thus 



Ehz — 



(AJ 



ab) 



(Nab)/ 2 



> C N/2 /(N/2), 



(IV.9) 



The experimental observation of Ehz < C n „/2/{no/2) 
implies that the two-mode entangled state consists of 
more than no atoms in total, that is, that there are nec- 
essarily more than n atoms entangled e.g. for n = 100, 
this could be a state like {|50)|51) -I- |51}|50)}/V2. As 
explained in the Appendix A, the result follows because 
values of Cj / J for increasing J are a decreasing sequence, 
and thus a measured value of Cj/J = 0.15, for example, 
can only come from an entangled state with at least 100 
atoms. 

The best HZ inter-mode entanglement is achieved in 
the attractive regime We note that the strongest theo- 
retical entropic entanglement e(p) j(32l |63| is found when 
all atom numbers are equally represented in the super- 
position. It is shown in [29] that the closest state to this 
optimum is obtained at a critical value of Ngu/n ~ —2.0, 
that is, the attractive interaction regime (as found in K 
and 7 Li isotopes) gives rise to a maximal spread in the 
distribution of numbers in each well. 

Interestingly, Fig. [7J shows the same point of maxi- 
mum is observed for the higher order entanglement mea- 

(2) 

sure E HZ . This measure can only detect entanglement 
that originates from superpositions of the type 1 50) 1 51) + 
1 5 1) 1 50) + 1 52) 1 49) , where at least some of the states 
of the superposition are separated by 2 quanta (Ap- 
pendix B). Similarly, the third order entanglement cri- 

(3) 

terion E HZ detects entanglement originating from states 
separated by 3 quanta. In the case of Fig [JJ where there 
is N = 100 quanta, the existence of entangled states such 
as |0> 1 100) + ... + 1100)|0) could be detected by measuring 



E^hz^ < 1- Higher order entanglement (e.g. E^ 1 ' < 1 
would not be possible. 



(101) 



3. Repulsive interactions 

The repulsive regime of positive g also predicts consid- 
erable planar squeezing (Fig. [Jj, but, in that case, the 
best planar squeezing is rotated into the X — Z plane as 
graphed in Fig. [5] [U, . A depiction of the resulting 
planar squeezing ellipsoid is shown in Fig 1101 

Thus, the corresponding HZ entanglement is between 
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Figure 9. (Color online) The repulsive interaction case for 
iV = 100, showing individual spin variances, and mean spin, 
{Jab)) f° r the ground state solution of Hamiltonian in the 
regime where there is a strong repulsive self-interaction g/n, 
for N — 100. Here k > is fixed and g is varied. 

the modes defined by the rotated coordinates, 

a' = (a + b) I (V2i) , b' = (a-b) /y/2. (IV.10) 




Figure 10. (Color online) The 3-D variance ellipsoid corre- 
sponding to N = 100 and a repulsive interaction at the opti- 
mum coupling of Nq/k = 40. Spin variances are reduced in 
both axes parallel to the X — Z plane, to show strong, but not 
perfect, planar quantum squeezing. The variance increases 
perpendicular to the squeezing plane, along the Y axis. 
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The corresponding entanglement criterion is given by: 



< E't 



HZ 



(AJ 



ab) 



<AJ Z ab) 



((ato) + (&t&))/2 



(IV.ll) 



The detection of spatial HZ entanglement between the 
two wells in the repulsive case would therefore require a 
different detection scheme, as proposed in [6l|. We note 
that in both repulsive and attractive cases, the HZ entan- 
glement can be very significant, so that the EPR steering 
nonlocality Eq (|III.13|) is predicted via measurement of 
both the first and second order HZ moments. Figure [§] 
indicates that, for fixed N, the repulsive case shows an 
increasing and then reducing first order HZ entanglement 
fill .f Op . as the nonlinearity g/n increases. The optimum 
case occurs for N = 100 and a repulsive interaction at a 
coupling of Ng/ k = 40. The squeezing ellipsoid for this 
coupling is shown in Fig [TO] 



use the above strategies to confirm an EPR-type entan- 
glement, one would measure the local EPR observables, 
Xa/b an d Pa/ Bi a t each well. This is because the mo- 
ments of (jIII.6[) are in terms of operators, a and b, which 
are linear combinations of the hermitian observables, X's 
and P's. Optically, the X and P are field amplitudes, 
the measurement of which is normally achieved by phase 
sensitive local oscillators [6fl. 



V. EPR ENTANGLEMENT: FOUR 
COMPONENT CASE 

We will show in this section how to use two additional 
modes per site to perform an effective "local oscillator" 
measurement in this BEC case. Similar strategies have 
been suggested by Ferris et al [l6| among others. 
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Figure 11. (Color online) Higher order entanglement for the 
case of N — 100. Other parameters as in [7] The red dashed 
line shows a reduced Ehz entanglement as g/n is increased 
above a certain level. The solid black line (starred) shows 

(2) 

the second order entanglement. Here E HZ < 1, indicating 
entanglement, is possible at higher g/n ratios where the first 
order (dashed) criterion shows no entanglement. 



Interestingly, however, from Fig. 111! we see that the 
second order entanglement criterion for TV = 100 picks 
up more entanglement, suggestive that the drop is due 
to a change in the nature of the entanglement, rather 
than to a loss of entanglement itself. Fig. [5] shows that 
a similar behavior occurs at much lower particle number 
with N = 6, although with less overall entanglement at 
the optimum coupling. 

The entanglement itself, via the measure of Ehz ■, can 
be deduced via measurements of the combined spins Jab , 
using interference measurements between the two con- 
densates, as has been performed in (24|. Results ob- 
tained in this fashion are important in confirming the 
existence of that type of entanglement within quantum 
theory, but as the measurements are not localized at each 
site, they cannot be viewed as rigorous tests of EPR en- 
tanglement, EPR steering or nonlocality. In order to 



A. Spin-operator entanglement criteria 

A true EPR experiment would involve coherent com- 
bination of second fields or condensates at each site, as 
depicted schematically in Fig. [5J To observe true EPR 
entanglement between sites A, B, a useful procedure is 
to use two modes per EPR site. Local intra-well spin 
measurements are defined: for well A. 



J A 



a\a 2 + a\aA /2, 



J A = [a[a 2 - a 2 ai ) /(2i), 



a\a 2 - a\aA /2, 



N A 



a\a 2 + a\a\. 



(V.l) 



Here ai i2 are mode operators for different components 
of the same site, typically different spatial modes or dif- 
ferent nuclear spins at each site. We will also introduce 
the notation for the corresponding raising and lowering 
spin operators, J A = J* + iJ\. Similar spin opera- 
tors are defined for site B. This defines complementary 
observables that are locally measurable at each site, us- 
ing Rabi rotations and number-difference measurements. 
Calculations of spin correlations at two sites can be car- 
ried out most simply on imaging on a micron scale, then 
dividing the imaged atoms into two halves for measure- 
ment purposes. A more sophisticated method is to add a 
time-dependent external potential to divide the conden- 
sate into two widely separated parts. While this gives 
results that depend on the potential, it provides a phys- 
ical separation between the sites. 

Having defined local spin operators, we now need to 
consider a suitable EPR entanglement measure. We 
present HZ-type criteria that are exp ressed in terms of 
these effective local spin operators [65l - l67| . Entanglement 
is confirmed if 



(J A Jb 



-\ |2 



> {JaJaJb^b^ 



(V.2) 
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This inequality thus uses operators which are measurable 
locally using Rabi rotations and number measurements 
[25| . Criteria involving higher moments are also possi- 
ble, but are not examined here. As for the original HZ 
criterion, the spin criterion can be rewritten using the 
procedure outlined in (29|. If we define Z = j\ J B , then 
we can easily show that A 2 (J^J^) = (Ja^a^b^b) ~ 

\(JlJ B )\ 2 - (JaJ a JbJb)<({Ja^1}JbJb) 

=2(Jp+J B ). (V.3) 
Similarly, defining Z' = J^Jg, one can show that 



{J\Jb)\ 2 >0. Thus, 



(J A Jb) 



(V.4) 



The spin entanglement criterion (|V.2[) becomes 

A 2 (JPb) 



E 



spin (1) 
HZ 



min{2(Jp+J^),2{J+JjJi)] < 1 ' (V - 5) 
i.e. HZ-type spin entanglement can then be verified if 
0<£f™ (1) < 1. 



Spin versions of the EPR paradox criteria piI.4IIIL5]) 
are possible (30l . IfM [6^ | , but here we focus on the spin 
EPR steering inequalities that relate to (|V.2j) . These 
have been derived in (70|: EPR steering is detected if 



\(JPb)\ 2 >([{Ja) 2 -{J Z a) 2 ±J Z a\ 



A(JbY-{JZY\) 



which can be rewritten as 



< E 



spin (1) 
HZ 



{Ja^a^b^b) \(Ja^b)\ 2 
rmn[2(Jp+J B ),2(J+J A J§) 



(V.6) 



1 

(V.7) 



We note the spin moments of Eqs (|V.5j) and (|V.7j) are 
actually measured via the x and y spin components, for 
example, using the expansion: 



{Ja^b) —{JaJb ~ iJaJb + J\Jb + J a Jb) 



(V.8) 



B. Linear multimode case 



We examine the linear case first, to model a fixed 
number of atoms with a minimal BEC nonlinear self 
interaction. Suppose a Fock number state \ipi n ) = 
l-^i)oi„i|^V2)o 4 „2|0)i>ini|0)6 <n2 is incident on a beam split- 
ter (Fig. [T2"|) , so that Ni and N% are fixed, and modes 
within each pair ai, b\ and a 2 , b 2 are coupled by the BS 
interaction, with a\ and 02 (and &i, b 2 ) remaining un- 
coupled. Output modes a\ and b\ are number-conserved 
according to ()II.2|) ; as is pair 02, b 2 , and are given as 

01,2 = (Oinl, 2+binl,2) / V2, and 6i i2 = (ainl,2-bi n x,2) / V2. 

The output state is 



iVi JV 2 



|out) = c n , n ,\n) a i\ri) a2 \Ni - n) b i\N 2 - n') b2 

n=0 n'=0 

(V.9) 



|N i: >|N 2 > 



|0>|0> 




|n> ai |n'^ 2 |Nx-n> bl |IW: 



b2 



Figure 12. (Color online) We consider pairs of Fock states 
transmitted through a beam splitter. Pair ai, fei are coupled 
and become entangled, as do <Z2, 62- 



where c n , n > = ^N^.N^./ y / 2 N ^+ N --n\{N 1 - n)\n'\{N 2 - n')\. 
We can evaluate moments, to obtain the prediction for 
the HZ spin criterion Eq. (|V.5j) . Fig. [T3l shows the result 
of varying N\ for fixed N% = 100. The asymmetric case 
is favorable to detecting entanglement. 
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Figure 13. (Color online) The entanglement of pairs of Fock 
states transmitted through a beam splitter can be detected 
via criterion Eq. (|V.5|) for the asymmetric case where pair 
a 2 and 62 have much greater numbers N2 3> Ni {Ni = 100). 
Entanglement is confirmed if S^f™' 1 ' < 1. 

Where the initial state is more complex, such as 
\ip in ) = lAMaiJA^kJ^ikJ^k^, the output 
state will involve superpositions of only even numbers 
of atoms in the symmetric and antisymmetric modes, so 

that I {J\J B ) l 2 H (a 2 aib 2 b{) | 2 = 0- As in t]i e case of 
Section IV. B. 2, we can detect this entanglement using an 
appropriate second order spin criterion. 



C. Nonlinear four component BEC case 

We now consider the EPR entanglement that can be 
generated and measured when the modes interact to form 
the four-mode BEC ground state. We focus on set-ups 
that will enable the four mode case to produce an EPR 
entanglement that is the replica of the two-mode HZ en- 
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tanglement, as displayed in Figures ([T UTTj) . In this case, 
the second mode at each site may be thought of as part 
of a measurement system. 



1. Four-mode BEC Hamiltonian 

We assume the two-well, four- mode system of Fig [5] is 
described by the Hamiltonian: 



E 



gijajalajai 



(V.10) 

We consider two modes at each EPR site A and B, with 
four modes in total, as shown schematically in Fig. [5] 
This corresponds to the two component per well exper- 
iments of [251 ]. and somewhat less closely to the multi- 
mode interferometry experiments of j27|. Depending on 
the exact configuration, the local modes at each EPR site 
can be independent (in which case local cross couplings 
gij are zero (gi2 = 0)), or not independent, as would be 
the case where the modes are coupled by the BEC self 
interaction term, so the couplings cannot be "turned off", 
as in the set-up of [25| . The coupling constant is propor- 
tional to the three-dimensional S-wave scattering length, 
so that oc , as in the two-mode case. For example, 
a typical value of the S-wave scattering length for 87 Kb 
is an = 100. 4ao, where ao is a Bohr radius. Zero cross 
couplings are likely require spatial separation of the two 
local modes, as might be achievable with four wells. 

The Hamiltonian (|V.10[) with k = K\ = K2 is based on 
the assumption that the second pair of modes ct2, 62 are 
coupled between the wells in the same way as the first 
pair <zi,6i, which implies similar diffusion across wells. 
The case where K2 = 0, K\ 7^ is possible where dif- 
fusion across the wells can be controlled, as where the 
local modes represent separate wells. We will examine 
the predictions for both cases. 



2. Symmetric tunneling case 

The BEC nonlinearity can enhance the entanglement. 
This is evident on comparing with the case of zero atom- 
atom interaction {g^ — 0), which corresponds to the re- 
sult of the linear beamsplitter model (Fig. [T2")) , and is 
indicated by the large red circles in the Figures [TilfTBI 
First, we examine the case of symmetric inter- well tun- 
neling with k — K\ — so there is complete symmetry 
between the nonlocal setups, but a variable local cross 
coupling gi2- Figure [Til shows entanglement using the 
HZ spin criterion Eq. (|V.5[) . for the ground state, for 
cases of both zero and strong local couplings gi2- Asym- 
metric atom numbers with TVi <C 7V 2 are required for the 
best entanglement, however, as shown in the inset of Fig. 
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Figure 14. (Color online) Entanglement using adiabatic cool- 
ing to ground state in a two- well potential, at T — OK, for the 
four mode model of Fig. [2] with a variety of local cross cou- 
plings gtj. Here k — Ki — K2 > and gn is varied with the 
other values of gij held in a fixed ratio, with Ni =5, N2 — 
100. E^z^ 1 ' < 1 indicates entanglement; E 3 ^^ < 0.5 in- 
dicates EPR steering. Curves are labelled in order of nesting 
as: (purple dash dotted) equal couplings; (blue dots) non- 
zero cross-couplings corresponding to 87 Kb Feshbach reso- 
nance with an = 100. 4ao, ai2 = 80.8ao, 022 = 95.5ao; (black 
dashed) without cross-correlations gi2 = 0; 322 = gn; (green 
solid curve) negative relative cross-coupling g\\.g\2 < 0. The 
inset shows the effect of increasingly symmetric atom num- 
bers. 



We note from Fig. [TJ] that the entanglement is im- 
proved by using a "local oscillator"-type approach, in 
which the second modes 02, 62 are independent of the 
first at each location (512 = 0) (being only combined 
at the spin measurement stage (jV.ll) ) and are of much 
greater numbers (iV 2 3> N\). In addition however, we 
note from the black dashed curve of Fig. [TS] that better 
entanglement is obtained if the second "local oscillator" 
pair <X2, 62 are also entangled optimally, as given by the 
critical point of the plots in Fig. [7] Thus, the optimal 



E y H ' z is at ^2522/^2 ~ —2.03 for the modes 02 and 62, 
and at N±gn/ ' k\ rs —2.1 for modes a\ and b\ (as shown in 
the inset of Fig. [T5]). The choice ^2322 ~ Nigu therefore 
gives enhanced EPR spin entanglement (red solid curve 
of Fig. [Til) . 

The minimum of E^^ 1 corresponds to the minimum 

achievable for the HZ entanglement E^ z ; this minimum 
is presented for the case N\ = 100 in Fig. [7] Better en- 
tanglement is thus achieved by increasing the number of 
atoms N±, provided the other constraints, that N2 ^> Ni 
and gn and 322 correspond to the critical choice for each 
mode pair, are satisfied, as shown in Fig. 1151 Analytical 
details are given in the Appendix. 

It is interesting that the case of approximately equal 
couplings gn — g 2 2 — 512 is generally less favorable for 
the HZ spin entanglement (Figure [T?]). This can be un- 
derstood if we rewrite the Hamiltonian (jV.10[) in terms of 
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Figure 15. Entanglement using adiabatic cooling to ground 
state in a two- well potential, at T = OK. Here k = ki = Kz, 
312 = 0, and both gw and 322 are varied so that Nign/ni = 
Af2g22/«2 • E^ z n ^ < 1 indicates entanglement; Ej!^ < 
0.5 indicates EPR steering. Main graph: (black dashed curve) 
Ni =5, N 2 = 100; (blue dotted curve) Ni = 20, N 2 = 100. 
The curves are for values of local coupling that optimize E^ z 
for each mode pair, in which case for N2 3> Ni the E a ^P z 
becomes the E^ z displayed in Fig. [Jj The inset reveals the 
individual degree of HZ entanglement E^ z for the mode pairs 
ai,6i and 02,62, as explained in the text. 



the spin operators. We obtain H ~ \ {JaY + X {JbY + 
K{a\bi+aib\+a\b 2 + a 2 b\), where \ — 5(511+922-2312) 
gives the effective nonlinearity, and those terms related 
to Jjf B , Nf 2 , N\ t 2 have been omitted. For equal cou- 
plings gi2 = pii = .922; the Hamiltonian thus effectively 
reduces to the linear term of the BS model of Fig. [T^l 
This is evident in the results of Figs. [H] and [TS] Fur- 
thermore, enhancement of the nonlinearity is possible, if 
g±2 becomes negative. The green solid curve of Figure 
[T4l shows an enhanced entanglement for negative local 
cross-coupling, gi2 < 0. 

The spin HZ entanglement is optimal in the attractive 
regime, gn < 0. Enhancement of entanglement in the 
repulsive regime is possible (Figure I16[) , if one examines 
the spin HZ entanglement for the rotated modes, a' , b' 
ofEq. PVTO)) . 

The effect of temperature is given in Fig. 1171 In our 
calculations, we account for effects of finite temperatures 
by assuming a canonical ensemble of p — exp[— if/feeT], 
with an inter- well coupling of ft/feg = 50nK. The critical 
temperature for EPR entanglement signature is shown in 

Fig. [13 

The predictions in this paper are based on the assump- 
tion that the total number N of atoms is fixed. Entangle- 
ment (Ejj Z = 0.5), though not EPR-steering, is obtain- 
able between the output ports of a beam splitter with 
a number (Fock) state input, in the absence of nonlin- 
ear coupling terms, as was shown in Sections IIV1 B and 
IVlB. However, for coherent state inputs, which have a 
Poissonian number distribution, this entanglement is not 
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Figure 16. (Color online) Adiabatic cooling to ground state 
in a two- well potential, at T — OK. Here n = ki = K2. Pa- 
rameters are as for Fig. 1141 but the entanglement parameter 
is calculated for the rotated modes a', b' of Eq. (|IV.10|l . For 
large N2, the strength of entanglement measure is enough to 
confirm EPR steering via the criterion Eq. (|V.7|l . 



possible [7ll |. and we draw the conclusion that number 
fluctuations will have an important effect on the entan- 
glement. We have studied this effect in Refs. [2^,|6l| and 
suggested in [6lj how one can adapt the HZ-entanglement 
criteria to reduce the sensitivity to atom number fluctua- 
tions. The effect of particle fluctuations on entanglement 
and precision measurement has been studied recently by 
Hyllus et al (72j. However, such modifications have not 
been studied in this case for the EPR steering type of 
nonlocality. 



3. Asymmetric tunneling case 

An alternative strategy more closely aligned to those 
used in optics is to consider K2 = , ti\ 7^ 0. In this 
case, the modes 02 and 62 are uncoupled and indepen- 
dent. If they are prepared in coherent states |«2) |/?2) 
(we take 0:2 = P% = a, where a is real), with a large, 

approaches the value given in 



the entanglement E^ 1 ^ 1 



the two-mode case, by E^' z . We explain this as follows. 
For independent modes, as shown by equation (|VI.7|) of 
the Appendix, the HZ spin entanglement criterion (|V.2[) 
becomes, upon assuming coherent states for 02 and 62, 



\(a[bi)\ 



2 a 4 > 



\a[a 1 b\b 1 ){l + a 2 f 



(V.ll) 



which we see will approach the required two-mode en- 
tanglement level in the limit of large a. Figure IT51 plots 
the result with finite numbers of atoms for the case of 
optimal E^ z which occurs at Nign/iti s» —2.03 when 
Ni = 100. We can see that the four mode EPR entan- 
glement achieved {Cj/J s» 0.15) is that of the two-mode 
case (Fig. [7]) provided there is a large enough number of 
atoms in the second mode. 
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Figure 17. (Color online) The effect of critical temperatures 
corresponding to the parameters of Fig. 1141 when Ng/n ~ 
-2.23. 
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Figure 18. (Color online) The effect of an uncorrelated co- 
herent atomic oscillator field for mode 2 in a coherent state 
with amplitude a, in the optimal case of Fig. [7] when 
Njg/K ~ -2.03. 



VI. SUMMARY 

We have examined a number of strategies capa- 
ble of generating detectable entanglement between two 
spatially-separated potential wells in a BEC. These in- 
clude both two and four-mode strategies similar to those 
already used for spin-squeezing, but generalized to a dou- 
ble well. The model used to calculate the relevant vari- 
ances has been shown to give a good fit to experimental 
data [25l . [fj3 | . Our results find that local cross couplings 
can have a strong effect on entanglement, and results for 
the EPR entanglement improve with higher atom num- 
bers. We find that a spin version of the Hillery-Zubairy 
(HZ) entanglement criterion appears readily suited to an- 
alyzing entanglement and the EPR paradox in these ex- 
periments. Furthermore, we have shown that the HZ en- 
tanglement criteria give information about the number 
of particles involved in the entangled state and, through 



higher order moments, the nature of the multiparticle 
entanglement. 
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APPENDIX A 

The system generally is described by a density matrix 
p = P R p R . Each p R is a pure state with a mean 
total number N R of atoms in modes a and b. The sums 
of the variances must be greater than the average of the 
variances of each component of the mixture [73j . The 
spin operators are defined according to (|IV.1[) . We now 
suppose the maximum number of atoms in any compo- 
nent to be Hq. Then the maximum spin for any such 
component is Jo = no/2. Then we are able to conclude 

that {{A R J^ B f + {A R j\ B ) 2 }/(N R /2) > CjJJo be- 
cause the functions Cj /Jq decrease monotonically with 
increasing Jq, as shown in the Ref. [59j . This leads to: 



(Aj* B y + (Ajx B y >Y,p R {(A R j$ B y + (A^n 

R 

^P R (N R /2)CjJJ (VI.1) 



> 



R 



Factorizing and using (TV) = J^r PrN r gives the final 
result 



r (AJf B ) 2 + (AJl B ) 2 

OffZ = — — > L.no/2/( n 0/ 2). 



{Nab) 12 



(VI.2) 



The observation of Ehz < C„ / 2 /(no/2) thus implies 
that the two mode entangled state consists of more than 
no atoms in total, that is, that there are necessarily more 
than no atoms entangled. In this case, the spin J AB 

measurement must be compatible with J^ B as defined 
in terms of just two modes, a and b. 



APPENDIX B 

The system generally is described by a density matrix 

P = Efl p RP?e P + Efl' p R'Pent> where pfep and p^ nt rep- 
resent pure separable and entangled states respectively. 
The higher-order HZ entanglement measure can therefore 
be written as a ratio 



R 



((at 



(6t) m b" 



(VI.3) 
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where 

R W 

and 

((at) m a m (fet)" 1 6™) = J2 Pr(^T « m faT & m >* 

R 

+ J2Pn>((*T* m (bTb m )R' 

R' 

Here (O) r represents the expectation value of O for state 
Pr. Since for a separable state, R < 1, we can see that 
if J2r> PR'( am (b ) )r.' = 0, it is always the case that p 
predicts R < 1. In short, the higher order entanglement, 
jr;( m ) < ]_ ( cannot be achieved unless there is a nonzero 
probability Pr< for a pure entangled state p^ nt for which 
(a m (6 T ) ) 7^ 0. Expanding the pure entangled state in 
terms of the number states yields most generally p = 
\ijj)(tjj\, where 

\<P) =J2cni\n)\l). (VIA) 

n J 

If only adjacent number states of type \n), \n + 1) have 
nonzero amplitude c„/, the (a m (6^) m ) = where m > 
1. Hence, the superposition (|VI.4[) necessarily includes 
states separated by m. 

APPENDIX C 

In this Appendix, we show we show how to directly 
"convert" the inter- well entanglement shown in Fig. [7] to 
an EPR entanglement, with the use of a "local oscillator"- 
type treatment which applies where two of the strong 
local modes are uncorrelated. This is the case of <?i 2 = 0, 
illustrated in Fig [2j 

Local oscillator measurements are achieved optically 
by combining a mode with a very strong coherent state 
[0| . We can achieve something effectively equivalent to 
a "local oscillator" measurement, where the second pair 



of levels a 2 , 6 2 are much more heavily populated than 
levels a\ and b\, by assuming the second pair of modes 
are in an uncorrelated coherent state. We explain this as 
follows. Since j\ = aja 2 and = a\a\ and Jg = &{6 2 
and Jg — bib\ we can rewrite the criterion (|V.2[) in terms 
of the mode operator moments, for this special case, by 
the factorization that is justified for independent fields 
at each location. Thus, 

\(JiJB)\ 2 = \(4bi)(a 2 bt)\\ (VI.5) 

and similarly 

((JaJaKJPb)) = (a\a 2 a ia lb\b2bibl) . (VI.6) 
The criterion (|V.2[) becomes 
|(e4&i)| 2 |(a 2 62)| 2 > ( a I a i 

6^ 1 )<(l + 4a 2 )(l + 6^ 2 )) 

(VI.7) 

Clearly, since the inter-well entanglement studied in 
Section HVl and summarized in Fig. [7]enables |(a|fei)| 2 > 
(a\aib\bi) via the HZ entanglement criterion, we will 
have (at least) the same level of four mode EPR entan- 
glement, provided 

|(a 2 4)| 2 >((l + 4a 2 )(l + 46 2 )). (VI.8) 

In fact, the inequality would represent violation of the 
two-site version of the Bell inequality discussed in [57}, 
which is not achievable for this system. However, it is still 
possible to optimize the EPR entanglement. This can be 
achieved in the following way. If the two modes a 2 and 
&2 are also coupled via an inter-well interaction (k 2 7^ 
in Fig. [2|), to produce the ground state solution of Fig 
[21 then Eflz < 1 amounts to Ka^)! 2 > (a^a^^)- 
The optimal E^ z is at A 2 (? 22 /k 2 —2.03, while for 
the modes a\ and a 2 , the optimal (|VI.7[) occurs for 
Nign/ni w —2.1 (inset of Fig. [IT)]) . This choice gives 
enhanced EPR entanglement as shown in Fig. [TJ] Bet- 
ter entanglement is possible for this optimal choice, as 
the numbers are increased (Fig. [15]) . 
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